Heat  Conduction  by  Molecular  Dynamics  Technique 


Sebastian  Volz 


National  Engineering  School  of  Mechanics  and  Aerotechnics 
Laboratory  of  Thermal  Studies  UMR  CNRS  6608 

Poitiers,  France 


Denis  Lemonnier  -  Lab.  of  Thermal  Studies  -  Poitiers 
Jean-Bernard  Saulnier  -  Lab.  of  Thermal  Studies  -  Poitiers 
Gang  Chen  -  NanoHeat  Transfer  and  Thermoelectrics  Lab.  -  UCLA 
Pierre  Beauchamp  -  Laboratoire  de  Metallurgie  Physique  -  Poitiers 


[my  A 


REPORT  DOCUMENTATION  PAGE  Form  Approved  0MB  No. 

0704-0188 

Public  reporting  burder  for  this  collection  of  information  is  estibated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing 
and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burder  to  Department  of  Defense,  Washington 
Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of 
law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  0MB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  EORM  TO  THE  ABOVE  ADDRESS. 

1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (EROM  -  TO) 

30-05-2001  Workshop  Presentations  30-05-2001  to  01-06-2001 

4.  TITLE  AND  SUBTITLE 

Heat  Conduction  by  Molecular  Dynamics  Technique 

Unclassified 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Volz,  Sebastian  ; 

Lemonnier,  Denis  ; 

Saulnier,  Jean-Bemard  ; 

Chen,  Gang  ; 

Beauchamp,  Pierre  ; 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PEREORMING  ORGANIZATION  NAME  AND  ADDRESS 

National  Engineering  School  of  Mechanics  and  Aerotechnics 

Laboratory  of  Thermal  Studies 

UMR  CNRS  6608 

Poitiers,  Erancexxxxx 

8.  PEREORMING  ORGANIZATION  REPORT 
NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME  AND  ADDRESS 

Office  of  Naval  Research  International  Eield  Office 

Office  of  Naval  Research 

Washington,  DCxxxxx 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 

NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

APUBLIC  RELEASE 

13.  SUPPLEMENTARY  NOTES 

See  Also  ADM001348,  Thermal  Materials  Workshop  2001,  held  in  Cambridge,  UK  on  May  30-June  1,  2001.  Additional  papers  can  be 
downloaded  from:  http://www-mech.eng.cam.ac.uk/onr/ 

14.  ABSTRACT 

UNDERSTANDING  AND  MONITORING  MATERIALS  PROPERTIES 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 
Public  Release 

18. 

NUMBER 
OF  PAGES 
21 

19.  NAME  OF  RESPONSIBLE  PERSON 

Fenster,  Lynn 
lfenster@dtic.mil 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 

Unclassified  Unclassified  Unclassified 

19b.  TELEPHONE  NUMBER 

International  Area  Code 

Area  Code  Telephone  Number 

703767-9007 

DSN 

427-9007 

Standard  Form  298  (Rev.  8-98) 
Prescribed  by  ANSI  Std  Z39.18 


UNDERSTANDING  AND  MONITORING  MATERIALS  PROPERTIES 


NEW  MATERIALS  contain  nano-micro  arch itectu red  structures 


o  NANOFIBERS: 

multiscale  complex  materials 
ultra  insulating  (0.007  W/mK) 


o  SUPERLATTICES 
monitoring  thermai  conductivity 
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o  NANOWIRES  templates: 


monitoring  anisotropy 


LOW-DIMENSIONAL  PHYSICS  FOR  HEAT  CONDUCTION 


o  Ballistic  Transport  of  Phonons 
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Boundary  Scattering 
More  resistive 

Specular 


SITUATION  of  MD 


^ _ Small  Scales 


honon  Transport:  Ballistic  - ►  Diffusive 

Phonon  Particle:  Boltzmann  TE 

No  interference,  Interface  transfer 
Small  and  Large  scales 


Phonon  Wave  ■ 

No  p-p  scattering,  limited  by  wl 
Interface  transfer  easier 


Molecular  Dynamics  -  MD 
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Atomic  Period  0.1  ps  Reiaxation  time  lOOps 
Lattice  Constant  O.Snm  Mean  Free  Path  lOOnm 


Ciassical  Heat  Conduction 
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MOLECULAR  DYNAMICS  TECHNIQUE 


COMPUTE  ALL  ATOMIC  TRAJECTORIES 


•  2nd  NEWTON  LAW 


•  STILLINGER-WEBER  POTENTIAL 


Coordonnee  X  (Angstroem) 
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ADVANTAGES  OF  MD  TECHNIQUE 


o  Phonon  Scattering  is  Difficuit  to  Model: 

Phonon  Particle  Approach:  Relaxation  Time 
Phonon  Wave  Approach:  No  scattering. 

MD  PROVIDES  A  COMPLETE  DESCRIPTION 
OF  PHONON  SCATTERING 

EXAMPLE:  NANOWIRE 

o  Phonon  Transport  Approach  Assumes  Fully  Periodic  Lattices 
MD  ALLOWS  TO  INCLUDE  ATOMIC  DEFAULTS  and  STRAINS 
EXAMPLE:  SUPERLATTICE 


o  Non-Equilibrium  Short  Time  Heat  Conduction 
MD  DESCRIBE  HT  BEHAVIOUR  AT  GigaHTz  FREQUENCIES 
EXAMPLE:  IN  BULK  SI 
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HEAT  FLUX  by  MD 


Kinetic  Term 
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THERMAL  CONDUCTIVITY  by  MD 


SILICON  NANOWIRE  MD  MODEL 


PERIODIC  BOUNDARY 


RIGID  BOUNDARY 


BOLTZMANN  TRANSPORT  EQUATION 
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S.  Volz  and  G.  Chen,  Heat  and  Technology,  18,  37,  2000. 
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COMPARISON  BETWEEN  MD&BTE  RESULTS 


l^t 


Is  boundary  scattering  the  only  cause  forthermal  conductivity  reduct 


on? 


S.  Volz  and  G.  Chen,  Applied  Physics  Letters,  57,  2056,  1999. 
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PHONON  CONFINEMENT  EFFECT  ON  HEAT  CONDUCTION 


o  Ballistic  Transport  of  Phonons 


RPTE  -  THE  DISCRETE  ORDINATE  METHOD 


QgradL,=;t-,(L:-4) 
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o  Ss,  discretizing  in  48  directions  (Q,  w^)  for  finite  iength  wire 

‘  Confinement  effect 

RIGID  -PERFECTLY  REFLECTING  BOUNDARY 
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RIGID  -PERFECTLY  REFLECTING  BOUNDARY 
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PHONON  CONFINEMENT  vs  BOUNDARY  SCATTERING 


PHONON  CONFINEMENT:  50%  REDUCTION 
BOUNDARY  SCATTERING:  70%  REDUCTION 


S.  Volz  and  D.  Lemonnier,  Physics  of  Low-Dimensional  Structures,  5/6,  91,  2000. 


Si/Ge  SUPERLATTICE  MD  MODELING 
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DISPLACEMENT  (A) 


STRAIN  EFFECT  ON  SUPERLATTICE  STRUCTURE 


o  Starting  with  mean  lattice  constant 
o  Implementing  Conjugate  Gradient  Method 


[001]  ATOMIC  PLANE  NUMBER 


CROSS-PLANE  THERMAL  CONDUCTIVITY  (W/mK) 


SUPERLATTICE  THERMAL  CONDUCTIVITY 


LAYER  THICKNESS  (A) 

^  S.  Volz,  J.B  Saulnier,  G.  Chen,  P.  Beauchamp,  Microelectronics  Journal  31  (9-10)  815,  2000. 
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EFFECTIVE  THERMAL  CONDUCTIVITY  AT  GIGAHERTZ  FREQUENCIES 


o  Fluctuation  Dissipation  Theorem 

bg  ^ 


q 


dt. 


gradT 


o  co  ''  dependence  at  Giga  frequencies 
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CONCLUSION 


MOLECULAR  DYNAMICS  TECHNIQUE: 

o  COMPLETELY  DESCRIBES  PHONON  SCATTERING 
o  ALLOWS  THE  SIMULATION  OF  DEFAULTS/STRAINS 
o  GIVES  ACCESS  TO  NON-EQUILIBRIUM  REGIMES 

-  IS  VERY  HEAVY  IN  TERMS  OF  COMPUTATION  TIME 

-  RELIES  ON  THE  INTERACTION  POTENTIAL  VALIDITY 


-  DOES  NOT  INCLUDE  QUANTUM  EFFECTS 


Effective  Thermal  Conductivity  (A.U.) 


ULTRA  SHORT  TIME  HEAT  CONDUCTION 
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ro-1  LAW  FOR  Si  THERMAL  CONDUCTIVITY  AT  GIGAHTZ  FREQUENCIES 
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